Resonance effects on the crossover of bosonic to fermionic superfluidity 
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Feshbach scattering resonances are being utilized in atomic gases to explore the entire crossover re- 
gion from a Bose-Einstein Condensation (BEC) of composite bosons to a Bardeen-Cooper-Schrieffer 
(BCS) of Cooper pairs. Several theoretical descriptions of the crossover have been developed based 
on an assumption that the fermionic interactions are dependent only on the value of a single mi- 
croscopic parameter, the scattering length for the interaction of fermion particles. Such a picture 
is not universal, however, and is only applicable to describe a system with an energetically broad 
Feshbach resonance. In the more general case in which narrow Feshbach resonances are included in 
the discussion, one must consider how the energy dependence of the scattering phase shift affects the 
physical properties of the system. We develop a theoretical framework which allows for a tuning of 
the scattering phase shift and its energy dependence, whose parameters can be fixed from realistic 
scattering solutions of the atomic physics. We show that BCS-like nonlocal solutions may build up 
in conditions of resonance scattering, depending on the effective range of the interactions. 

PACS numbers: 03.75.Ss 



Introduction 

Fermion pairing is a fundamental concept in the man- 
ifestation of non-trivial ground states in condensed- 
matter physics. Condensation of composite fermions, 
viewed as bound states of an electron and an even num- 
ber of vortices, is a possible explanation of Integer and 
Fractional Quantum Hall Effects in the highly degenerate 
two-dimensional electron gas, the vortices corresponding 
to the fractionally charged quasi-particles of Laughlin's 
theory pj . Exciton formation and possibly Bose-Einstein 
condensation 0] in semiconductor structures is one more 
example Cooper pairing resulting from correlations 
in momentum space is the mechanism determining the 
BCS-type superconductivity in metallic compounds 0, 
while strong correlations [j| and real-space pairing with a 
short coherence length characterize the high-temperature 
superconductors, where pair correlations manifest in the 
opening of a pseudogap well before the superconducting 
transition 

The achievement of quantum degeneracy in atomic 
Fermi gases of 40 K and 6 Li after cooling in dipolar optical 
traps |3j la ISl lid - llll I12J ) and of accurate control of the 
interactions by means of Feshbach resonances ^3|j nas 
made these concepts accessible in atomic-physics experi- 
ments. This has initiated an exploration of an intriguing 
system which should reveal the predicted resonance su- 
perfluidity [Til ITU fill and the nature of BCS pairing 



with large attractive interaction [T?|. 

In these experiments, the magnetic field is changed to 
tune the position of a bound state in a closed channel 
in the interatomic potential with respect to the thresh- 
old of zero scattering energy. The zero energy is defined 
by the asymptotic value of an open channel potential in 
which the colliding atoms enter. At exact resonance, the 
bound state connects with the zero energy scattering so- 
lution and the scattering length is infinity. Variations 
of the magnetic field B with respect to the resonance 
value Bq can be converted into an effective detuning 
u=(B- B )A/x mas , with A^i mag the relative magnetic 
moment between the closed-channel bound state and the 
open-channel threshold. As displayed in Fig. 2] for the 
case of 6 Li, positive (negative) detunings with respect to 
the resonance correspond to effectively attractive (repul- 
sive) interactions resulting in negative (positive) scatter- 
ing lengths a. The Feshbach mechanism evidently in- 
volves a separation of energy and length scales between 
the background and the resonant behavior, the former 
driven by the value cib g at very large v, and the latter 
dictated by the position and width Ai/ = Ai?A/i mag of 
the resonance. We point out that certain physical sys- 
tems have large values for the background scattering due 
to the presence of a potential or shape resonance in the 
open channel, and a multiple resonance model should be 
developed in that case 

The side of the resonance corresponding to repulsive 
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interactions (a > 0) gives rise to a rich and complex 
quantum system. Here, the formation of weakly bound 
molecular states of two fermions has been obtained 
|2ll E3| and their Bose-Einstein condensation (BEC) [23, 
|24L 12a , l2q has been observed. On the other side of the 
Feshbach resonance (a < 0) where the interactions are 
attractive, experimental measurements are being made 
on the superfiuid paired state that is likely associated 
with a strong-coupling version of the Bardeen-Cooper- 
Schrieffer (BCS) theory. 
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FIG. 1: The broad Feshbach resonance in 6 Li [T^. The scat- 
tering length a is shown as a function of magnetic field, for 
the (f,m f ) = (1/2, -1/2) and (1/2, 1/2) mixed spin channel. 

By continuous tuning of the interactions across the Fes- 
hbach resonance, the whole crossover region from the rel- 
atively simple limits of a BEC of preformed pairs and 
a BCS-like superfiuid state is accessible in the exper- 
iment. There are many open questions in the under- 
standing of the region in between, including the relevant 
correlations in the system and the nature of the order 
parameter j^]. Several theoretical descriptions of the 
BCS-BEC crossover in Fermi gases with Feshbach reso- 
nances [23, Hi£ bM have been developed, inspired by 
the early seminal works by Leggett [12 and Nozieres and 
Schmitt-Rink |3^ |. and subsequent developments in the 
context of high-T c superconductivity. 

The crossover theories can be divided in two gen- 
eral classes: studies in which the interactions are 
parametrized by the scattering length [34| , and those that 
explicitly include the presence of molecular states of two 
fermions and the hybridization of these paired states with 
their unpaired counterpart [35j. Both approaches share 
the idea that the formation of Cooper pairs and their con- 
densation to the coherent superfiuid state do not occur 
at the same time. There are formal connections based 
on either application of the Hubbard-Stratonovich trans- 
formation to the functional integral [36| , or to the dress- 
ing of the open and closed channel solutions which 
demonstrate the equivalence of the approaches when the 



resonance state has a sufficiently short lifetime. 

Two relevant questions arise in the quantitative de- 
scription of the crossover behavior in atomic gases. First, 
an adequate many-body theory providing the crossover 
order parameter must be able to recover the results of 
the four-fermion scattering problem entering the Fesh- 
bach physics [3]} ■ Second, the energy dependence of the 
scattering phase shift may be important over the energy 
range of the Fermi energy [3^ , and thus the width of the 
resonance becomes a relevant parameter of the theory in 
addition to a. In the unitarity limit \a\ — > 00, where 
the scattering length is meaningless, the universal be- 
havior of the thermodynamic properties emerges only in 
the case of broad resonances 39] , and otherwise the spe- 
cific details of the resonance will enter. All these aspects 
pose stringent conditions to the choice of suited theoret- 
ical approaches. Along these lines, numerical methods 
would be a very useful guideline to assess the validity of 
approximate schemes and their underlying physics. 

In this work we focus on the issue of the energy de- 
pendence of the scattering phase shift in the unitarity 
limit, anticipating that the general case of a BCS-BEC 
crossover driven by two independent parameters (the de- 
tuning and width of the resonance) will lead to an in- 
triguing and nontrivial phase diagram. Along these lines, 
we develop a model that is able to interpolate between 
the two limits of a broad resonance, where the scatter- 
ing length resulting from low-energy resonance renormal- 
ization completely encapsulates the interaction behavior, 
and of a "high-quality" resonance in a broad Fermi sea, 
where large values of scattering energies need to be sam- 
pled. We model an interaction potential that is composed 
of a short-range attractive well followed by a barrier 
structure, where scattering length and resonance width 
can be independently tuned. We use such a model poten- 
tial to show that a nonlocal BCS-like superfiuid ground 
state may emerge under non trivial conditions. The re- 
sults are indicative at this stage, as they are obtained 
within a mean-field approach that is not necessarily valid 
in the regime of strong interactions. However, it is a nec- 
essary precursor to perform such calculations to guide the 
more complicated Quantum Monte Carlo (QMC) simu- 
lations, that are under way. 

The structure of the paper is as follows. We first de- 
velop an explicit microscopic model based on a potential- 
well and a potential-barrier to encapsulate the essen- 
tial physics of the Feshbach resonance. We then de- 
rive the nonlocal BCS equations for this system, assum- 
ing a closure of the hierarchy of many-body correlations 
at the BCS level. Finally, the BCS equations are self- 
consistently solved and the results discussed in view of 
their application to the development of a more detailed 
QMC simulation needed to capture the effects of many- 
particle correlations. 
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A well-barrier model of the Feshbach resonance 

We model the presence of a Feshbach resonance by the 
interaction potential of the form displayed in Fig. |3 
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V(r) = { Vx 




r < r o 

r < r < n 

otherwise 



(1) 



that is characterized by an attractive well with depth Vq 
and width ro and a barrier with height V\ and width r± — 
ro = r w . This well-barrier model allows us to incorporate 
the essential energy-dependence of the scattering physics 
into a single-channel scattering scheme. 
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different regimes of narrow and broad resonance based 
on a comparison of the Fermi energy Ep = h 2 kp/2m 
with kp — (S^n) 1 / 3 with the resonance width Av. The 
resonance width Av can be expressed in terms of the ma- 
trix element g for the coupling between the closed and 
open channels as Ai> — g^/n, that enters the resonance- 
superfluidity Hamiltonian, as pointed out in Ref. |4fj|. 

The parameters of the model potential can in princi- 
ple be adjusted to reproduce the scattering properties 
of an atomic sample, as e.g. they are determined from 
collision experiments or from full coupled-channel calcu- 
lations 0. We instead perform here a study over a set 
of parameters chosen in order to illustrate the important 
physical behavior of the system. We fix the scattering 
length to a large and positive value, a = 5000 ao, with 
ao the Bohr radius. We ensure this is large compared to 
the interparticle spacing as determined from the density 
which we take a,s n — 1.054 x 10 14 cm -3 . The large pos- 
itive scattering length corresponds to the region of uni- 
tarity limited behavior considered just on the BEC-side 
of the resonance. With these fixed constraints we vary 
the resonance width crossing the full region of broad to 
narrow values. To this aim, we first set the range ro of 
the potential to a value that is small enough to satisfy the 
condition nr^ -C 1. This is the case for r$ = 2000 ao giv- 
ing nr-Q = 0.125. We intentionally choose this value here 
so that it is not too small which will allow us later to see 
on a single energy scale both the scattering and many- 
body effects. Then, we solve for V\ (the barrier height) 
and n (the barrier width) which are not independent 
and can be tuned in order to change the tunneling rate 
through the barrier and therefore Av. 

This requires solving the scattering problem for V(r). 
The two-body scattering function ^(r) is given in the 
three regions by 



FIG. 2: A well-barrier model that allows for the independent 
tuning of the scattering length a and of the resonance width 
g. The well can support resonant states with a width de- 
pendent on the tunneling through the potential barrier. The 
situations shown correspond to the values used in subsequent 
simulations which all have a = 5000ao, where ao is the Bohr 
radius. Different cases are represented by the solid, dashed 
and dotted lines correspond to increasing values of the reso- 
nance width g (cases 1, 4, and 8 in Tabic 0. 

The condition for diluteness of the gas is that the range 
of the interatomic potential should be much less than the 
interparticle spacing, a condition that we assume to be 
fulfilled through this paper. Explicitly this requires us 
to satisfy at all times nr^ <C 1, where n is the particle 
density. It is important to emphasize that this is dis- 



tinct from the condition 



•C 1 associated with the 



convergence of the perturbation theory of the gas in the 
contact scattering approximation. In the vicinity of the 
resonance, we may have the situation in which |a| > ro so 
that we simultaneously have a system which is dilute in 
the sense that nrj) <C 1 but unitarity limited in the sense 
that no 3 > 1. We can furthermore distinguish the two 




r < r 

r < r < n 

otherwise 



(2) 



with 



r*/(r) 



sin(fc^r) 

A\ exp[-/c fc (r - r )] + B\ exp[fc b (r - r )] 
A-i sin(fcr) + B 2 cos(£x) (3) 



and 



k w = ^Jk 2 +mV /n 2 

h = 



\J mVi/ft 2 — k 2 



(4) 



The values of and B\,2 are determined after im- 
posing the usual boundary conditions ^ w (ro) = 'l'b(ro), 
*w(n>) = tf'tfo). *i(ri) =_*/(ri), and %(n) = 
ty'f(ri), i.e. ensuring the continuity of the wave func- 
tion and of its derivative at the boundaries ro and r\. 
The two-body T-matrix can then be extracted from the 
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scattering solution 



T(k) = 



4nh 2 B 2 



mk{iB 2 - A 2 ) ■ 



(5) 



We can now proceed to link the parameters of the 
potential with the parameters of the many-body theory 
we wish to solve. The first parameter, the scattering 
length a is determined simply from its definition applied 
to Eq. J5J, 



4nTi 2 a 



= lim T(k). 



(0) 



The width of the resonance is related physically to a ma- 
trix element between a continuum scattering state and a 
high-quality resonant state in the inner well. In practice, 
such a matrix element is determined by the tunneling 
rate through the barrier. In order to make the link ex- 
plicit we present the usual relation between the T and 
5-matrices in scattering theory 



T{k) 



2nfi 2 i 



[S(k) 1] 



(7) 



and give the Feshbach form for the S'-matrix in the one 
resonance parametrization 



S(k) 



2i%| 5 



4-Kh 2 r _ h 2 k 2 \ 
m ^ Tn ' 

Solving Eqs. Q and (JSJ) for g gives 



ik\g\ 



(8) 



m 2 R cB 



(9) 



where we have removed the background dbg = 
and the effective range is by definition R Q g = 
— (4irh 2 /m)(d 2 T(k)~ 1 /dk 2 )\ k=0 which is determined ex- 
plicitly from the scattering model Eq. JSJ . Note this 
form in Eq. © agrees with the expression in [j^. The 
dimensionless width in units of the Fermi energy is then 
the ratio g^/n/Ep as previously discussed. 

We list in the second column of Table |U the set of 
g^/n/Ep values that, together with the fixed condition 
a = +5000 ao, define from here on our case studies. The 
corresponding values of ri and Vq are reported in the 
third and fourth column. We illustrate in Fig. the 
wavefunctions r\E'(r) that are found for the three V(r) 
potentials reported in Fig. [5] and labeled as cases 1 , 8 
and 10 in Table [H and for the incoming kinetic energy 
Ek/ks = 1 nK. While it is expected that the wavefunc- 
tion should vanish near r — +5000 ao, corresponding to 
the value of a, it is clear that a knee develops in the 
probability amplitude for the scattered wave in the spa- 
tial region within the barrier. This feature becomes more 
pronounced as the resonance width decreases. 



TABLE I: Parameters for the model potential V(r) in Fig. [5] 
The values of r± and Vb are reported in the third and fourth 
column, that yield the different values of g^/n/EF in the sec- 
ond column with a — 5000ao fixed. The various cases are 
numbered as in the first column. The remaining model pa- 
rameters are ro = 2000ao and Vi = WOfiK. 



Case 


g^/n/Ep 


ri/ro 


V {pK) 


1 


0.94 


1.3345 


45.0581 


2 


0.99 


1.3236 


44.9236 


3 


1.05 


1.312 


44.7673 


4 


1.12 


1.2994 


44.5806 


5 


1.32 


1.2712 


44.0898 


6 


1.48 


1.255 


43.7547 


7 


1.71 


1.237 


43.3292 


8 


2.09 


1.217 


42.7804 


9 


2.96 


1.1935 


42.0159 



1.2 



Vn/E F = 0.94 
Vn/E F = 2.09 
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FIG. 3: The two-body wave function r^t(r) corresponding 
to the case potentials V(r) 1 and 8 in Fig. [5] and in Tab. [I] 
g^/n/E F = 0.94 (solid line), 2.09 (dotted line). The incoming 
kinetic energy is Ek/ks = 1 nK. 



The nonlocal BCS equations 

We can now utilize our analysis of the two-body physics 
for this model potential as an input to the many-body 
theory. In particular we wish to find the ground-state 
of the Fermi gas interacting via V(r). As mentioned in 
the introduction, we aim in this paper to solve the vari- 
ational BCS scheme for the case of a nonlocal potential 
interaction. Since our nonlocal potential allows us to ex- 
plore the consequences of scattering resonances, this is a 
suitable foundation as eventual input to QMC schemes. 
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Such schemes are necessary to include the many-particle 
correlations which go beyond the pairing fields studied 
here at the mean-field level in order to determine the 
unitarity limit properties of the gas. 

The ground state is determined within the variational 
BCS scheme by means of the wave- function 



|*o >= n k (w k + UkOk iT ol k) j_)|0 > 



(10) 



where 



is the creation operator for electrons of spin 
cr. The normalization of $o leads to the condition |u k | 2 + 
|wk| 2 = I- The expected value of the ground-state energy 
reads: 

E = 2e k |^k| 2 + Vk.k'Mk^k^k'^k' 

k kk' 

+ Vk.k'ttk^k^k'^k' en) 

kk' 

where ek = ?i 2 k 2 /2m. 

The summations in Eq. (|llf> can be converted into one- 
dimensional integrals after performing the angular inte- 
grations. For the generic function F(k'), in our situation 
of isotropic pairing this amounts to a substitution 



The second and third terms in Eq. (|16f) are the Hartree- 
Fock corrections to the single-particle self-energy, that 
manifest as a correction to the chemical potential. 

We iteratively solve Eqs. (|H|l and (|T5|l until self- 
consistency is obtained. Since two different length scales, 
ro and kp 1 , are involved, integrations on Eqs. I|14|) and 
(|15fl are executed by means of Gaussian integration in 
three different grids of fc-points that are in the ranges of 
fc-Fi ?-q~ j ano - (t ne latter is performed after changing 
to the inverse variable 1/fc). We have typically used a 
total of Nk = 1200 grid points for this purpose. 

The Bogoliubov Uk and Vk functions resulting from the 
energy minimization arc finally evaluated through the ex- 
pressions 




E(k) 



(17) 
(18) 



Self-consistent Results 



E^k'F(k') 



(2tt) 3 



dqq 2 V(k,q)F{q) , (12) 



with V(k, q) determined from the three-dimensional 
Fourier transform of the spatial potential 



27T 

V(k,q) = —(Vo + Vi) 



smr Q \k + q\ sinr |fc — q\ 



Vi 



sin n I k 



\k + q\ \k-q\ 
q\ sinrijA: — q\ 



1\ 



9 



(13) 



The BCS solution is obtained by minimizing the free 
energy / = ($ |#|<I>o) - (J,($ \N\$ ) with respect to 
the variational parameters Uk^v^. The chemical poten- 
tial /z is determined by the constraint to have the cor- 
rect particle density. The two resulting equations to 
be self-consistently solved correspond respectively to the 
isotropic superfiuid gap and the particle density and are 
given by 



A(fc) 



(2 



TT^J \ E(k)J 



(14) 
(15) 



In Eqs. l(T4^l - lfT3|) . the excitation energy E(k) = 
J A(/c) 2 + £(fc) 2 is expressed in terms of the gap func- 
tion and of the single-particle energy £(k) 



£(k) = £k - ^ + \v k =q=Qn - ^2 V k,v ( 1 



E(k<) 



(16) 



< 




FIG. 4: A(fc) in the unitarity limit with a = 5000 ao, as a 
function of the resonance width g^/n/EF- Curves with in- 
creasing values of A(fc = 0) correspond to decreasing values 
of gypnj Ef, from case 9 to case 1 in Table |H 

The gap function A(fc) is displayed in Fig. 0] for differ- 
ent values of g^fn/ Ef- The damped oscillatory behavior 
on the scale of 1 /rg is a manifestation of the pairing po- 
tential, with a wavelength and damping coefficient that 
are almost independent of the resonance width. 

The behavior of A = A(fc = 0) as a function of 
gy/n/Ep is summarized by the squares in Fig. [SI together 
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with the behavior of the chemical potential (circles) fol- 
lowing the gap. Fig. [5] demonstrates that a BCS-like so- 
lution emerges while the resonance shrinks on the scale 
of Ep, even on the BEC-side of the resonance, where the 
BCS- variational ansatz is not expected to give a complete 
description. These results suggest that a high-quality 
resonance leads to a nonuniversal regime in the unitarity 
limit. We have checked the consistency of this picture 
by computing the self-consistent solutions to Eqs. i|14|) - 
dl5fl in the case of narrower wells. The narrowing of the 




o 1 a 

gVn/E F 



FIG. 5: Ao = A(k — 0) (squares) and fj, (circles) as functions 
of g^Jn/Ep, with a = +5000 ao fixed. 

resonance seems to increase the level of the interactions 
that are responsible for the superfluid pairing. This is 
quantified in the ground-state energy whose values 
are seen to become larger and negative with decreasing 
g^/n/Ep (see Tab.|ll|l- 

TABLE II: Values of the total energy Eq (first column), chemi- 
cal potential fi (second column) and of the gap Ao = A(k — 0) 
at k = for the different values of g^/njEp in Tab.|l] 



g^/n/Ep 


Eo/Ep 


H/Ep 


Ao/kp 


0.94 


-1.286760 


2.757061 


1.967685 


0.99 


-1.240605 


2.661722 


1.838010 


1.05 


-1.191324 


2.563657 


1.699060 


1.12 


-1.138285 


2.461363 


1.545724 


1.32 


-1.021072 


2.247637 


1.195665 


1.48 


-0.953917 


2.133496 


0.993021 


1.71 


-0.877986 


2.012829 


0.771605 


2.09 


-0.816665 


1.910900 


0.556854 


2.96 


-0.687847 


1.751602 


0.330776 



Further insight can be obtained from the analysis of 
the momentum distribution n(k) 

n{k) = \v k \ 2 , (19) 

that is reported in Fig. In agreement with the conclu- 
sions from Fig. [5] larger values of g^/n/Ep correspond 
larger values of the jump at fcp and thus to a normal 
Fermi-gas character. The momentum distribution typi- 
cal of a BCS superfluid develops with decreasing values 
of gy/n/Ep. 

1 . 5 I — i i i i I i i i i I i i i i i i i i i 




0.5 1 1.5 2 

k/k F 

FIG. 6: Momentum distribution of the Fermi gas in the uni- 
tarity limit with a — 5000 ao, as a function of the resonance 
width g^/njEp. Curves with increasing values of n(k — 0) 
correspond to increasing values of g^/n/Ep, from case 1 to 
case 9 in Tab.Q] 

The superfluid state is also signaled by the emergence 
of a peak in the pair distribution function, that is 

g(r) =9HF(r)+g P (r) , (20) 
where the Hartree-Fock contribution is 

^w^-lW)^) 2 /* 6 ^' 1 "' 12 (21) 

and the pairing term 

5p(r) = + (si) / dke ~ 4k rM k«k . (22) 

The peak emerges on the scale of the interparticle dis- 
tance r s ciQ with r s — (47rnaQ/3) -1 / 3 , and becomes better 
defined as long as the quality of the resonant mode in- 
creases. 
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i 1 r 




0.5 1 1.5 2 

r/(r s a ) 

FIG. 7: Pair distribution function of the Fermi gas vs. the 
interparticle distance r s ao in the unitarity limit with a — 
5000 ao. Curves with increasing values of peak height corre- 
spond to decreasing values of gy/n/ Ef, from case 2 to case 8 
in Tab. □ 

Conclusions 

We have determined a solution to the many-body the- 
ory for the case of a Feshbach resonance within the BCS- 
variational ansatz. The situation we considered had a 
positive large scattering length corresponding to a near 



unitarity limited system on the BEC side of the reso- 
nance. The dependence of the scattering phase shift on 
energy was shown to be important in the self-consistent 
solutions for the mean- fields. The solutions were found 
to depend sensitively on the resonance width. We could 
vary all of these parameters within our model by mod- 
ifying the properties of an explicit and simple poten- 
tial model consisting of a potential well and potential 
barrier. Due to the fact that the solution we obtained 
was for a nonlocal system, the resulting many-body the- 
ory did not suffer from a formal ultraviolet divergence 
and was therefore automatically renormalized. The self- 
consistent pairing field had to be determined at each 
wavevector value, illuminating the nonlocal character of 
the solution. The fact that we find a BCS solution for a 
positive value of the scattering length can be understood 
from the total self-energy of the system, which can be 
negative for narrow Feshbach resonances. 

We emphasize that around the resonance, the BCS- 
variational ansatz is not expected to give a complete de- 
scription, although this is the exclusively applied frame- 
work in which the BCS-BEC crossover theories have been 
implemented to date. It is a necessary precursor to ex- 
plore and characterize the solutions we have found at the 
mean-field level as inputs to a Quantum Monte Carlo nu- 
merical method, which is able to encapsulate the many- 
particle correlations which we have dropped. This is an 
important problem where it is anticipated that the uni- 
tarity limit of the quantum system will emerge when the 
scattering resonance is broad, but that a high-quality res- 
onance will lead to a nonuniversal regime which depends 
on microscopic resonance parameters. 
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